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ABSTRACT 


The  onset  of  transition  from  laminar  to  turbulent  flow  in  a high  Reynolds  number, 
water  boundary  layer  has  been  investigated  with  emphasis  on  predicting  the  global 
effects  of  distributed  surface  roughness  on  transition.  The  particular  approach 
which  has  been  followed  is  to  use  the  resuits  of  linear  stability  theory  ~s  an 
indicator  of  transition,  but,  for  the  initial  results  which  are  reported  herein, 
only  comparisons  of  the  stability  results  themselves  (not  actual  transition  pre- 
dictions) are  given.  The  effects  of  surface  roughness  on  transition  have  been 
included  by  means  of  an  existing  phenomenological  model  for  the  effects  of  dis- 
tributed roughness  on  the  mean  flow  profiles.  This  model,  which  was  originally 
developed  for  high  Mach  number,  compressible  boundary  layers,  has  been  applied 
without  change  to  the  present  incompressible,  water  boundary-layer  environment, 
and  identical  values  of  the  required  empirical  constants  have  been  used  in  both 
cases.  Some  important  conclusions  which  can  be  obtained  from  the  numerical 
results  are  that  the  presence  of  favorable  pressure  gradients  and/or  surface 
heating  serve  to  make  the  boundary  layer  more  susceptible  to  roughness,  not  only 
when  compared  in  terms  of  the  actual  roughness  height,  k,  but  also  when  compared 
in  terms  of  the  roughness  height  to  momentum  thickness  ratio  (k/e).  Further  when 
the  roughness  is  sufficiently  large,  its  presence  can  change  the  effect  of  surface 
heat  addition  from  a strongly  stabilizing  factor  to  a strongly  destabl i 1 izing 
phenomena,  indicating  that  heating  on  walls  of  sufficiently  large  roughness  can 
be  detrimental  rather  than  helpful. 
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INTRODUCTION 


When  a high  Reynolds  number  vehicle  moves  through  a quiescent  fluid,  viscous 
forces  in  the  boundary  layer  accelerate  the  fluid  near  the  vehicle's  surface 
and  leave  it  with  a residual  net  momentum.  For  vehicles  on  which  the  boundary 
layer  remains  attached,  this  momentum  represents  nearly  the  entire  drag  of  the 
vehicle.  The  total  amount  of  momentum  which  is  transferred  to  the  surrounding 
fluid  (and,  hence,  the  drag  of  the  vehicle)  depends  crucially  on  the  character 
of  the  boundary  layer.  If  the  boundary  layer  is  laminar,  the  amount  of  momentum 
deposited  in  the  boundary  layer  is  relatively  small,  but,  if  the  boundary  layer 
is  turbulent,  it  is  substantially  larger.  Since  this  drag  must  be  overcome  by 
the  propulsive  system,  it  behooves  us  to  minimize,  to  as  large  an  extent  as 
practical,  the  amount  of  momentum  which  is  transferred  to  the  surrounding  fluid. 
One  particularly  attractive  method  for  achieving  this  objective  is  to  employ 
boundary  layer  control  techniques  to  delay  transition  from  the  laminar  to  the 
turbulent  state.  If  the  laminar  boundary  layer  can  be  maintained  over  the 
entire  vehicle  (or  nearly  so),  impressive  performance  improvements  can  be 
obtained. 

In  the  past  few  years,  substantial  progress  toward  reaching  this  goal  of  a high 
Reynolds  number,  all-laminar  vehicle  has  been  reported.  The  incorporation  of 
recent  advances  in  boundary  layer  control  techniques  has  allowed  the  feasibility 
of  all-laminar,  underwater  vehicles  to  be  demonstrated  at  practical  body  Reynolds 
numbers,  at  least  in  carefully  controlled  experimental  situations.  The  ability 
to  maintain  these  impressive  performance  improvements  in  a realistic  environment 
or  to  extend  them  to  even  higher  Reynolds  numbers  represents  an  important  achieve- 
ment which  must  still  be  demonstrated.  The  present  report  describes  some  initial 
results  of  a study  which  is  aimed  at  determining  how  vehicles  whose  design  is 
based  on  these  advanced  boundary- layer  control  techniques  will  perforin  in  such 
realistic  environments.  In  particular,  the  results  in  this  report  are  centered 
on  estimating  appropriate  levels  of  surface  roughness  which  can  be  tolerated 
without  negating  the  effects  of  the  boundary  layer  control. 

The  two  principal  boundary  layer  control  techniques  which  have  been  relied  upon 
to  delay  transition  in  the  experimental  demonstrations  noted  above  are  body 


shaping  and  surface  heating.  Body  shaping  is  used  to  generate  a favorable 
(negative)  pressure  gradient  over  large  portions  of  the  vehicle,  thus  providing 
a boundary  layer  which  is  strongly  stable  to  disturbances.  Surface  heating  is 
likewise  strongly  stabilizing  in  water  boundary  layers  because  of  the  simultane- 
ous effects  of  the  variation  of  the  viscosity  of  water  with  temperature  in  con- 
junction with  the  relatively  high  Prandtl  number  of  water.  Experimental  results 
have  consistently  shown  that  techniques  such  as  these  which  serve  to  increase 
the  stability  of  the  boundary  layer  also  serve  to  increase  the  transition  Reynolds 
number.  Some  computations  which  show  the  magnitude  of  these  two  factors  on 
stability  and  transition  have  been  given  by  Gazley  et  al.  (1976). 

Quite  obviously,  any  real  surface  is  not  perfectly  smooth  but  contains  some 
roughness.  It  is  well  known  that  surface  roughness  has+ens  boundary-layer  tran- 
sition, and,  as  will  be  shown  herein,  both  favorable  pressure  gradients  and  sur- 
face heating  tend  to  amplify  the  effects  of  roughness.  Surface  roughness  may 
enter  from  a variety  of  sources,  including  imperfections  arising  from  the  manu- 
facturing process  and  surface  degradation  with  time  due  to  external  influences. 

It  is  of  fundamental  importance  to  assess  the  interactions  between  surface  rough- 
ness and  body  shaping  and  heating  with  regard  to  their  effect  on  transition.  For 
example,  it  is  anticipated  that  there  should  be  a threshold  roughness  below  which 
the  transition  Reynolds  number  is  indistinguishable  from  that  on  smooth  walls. 

It  is  important  to  determine  whether  this  "smooth  wall"  range  is  wide  or  narrow, 
and  whether  the  effects  of  roughness  become  evident  gradually  or  catastrophically 
once  this  level  has  been  exceeded.  This  roughness  behavior  can  have  important 
practical  implications  as  regards  the  implementation  of  shaping  and  heating  in 
real  design  problems  where  the  degree  of  smoothness  which  is  practically  attain- 
able is  of  concern.  Further,  it  is  important  to  identify  acceptable  roughness 
heights  so  that  overly-conservati ve  surface  finishes  are  not  specified. 

The  particular  emphasis  of  the  present  report  is  to  make  a qualitative  assessment 
of  the  mutual  interactions  between  the  effects  of  surface  roughness  and  favorable 
pressure  gradient  and  surface  heating.  The  basic  approach  which  is  followed 
herein  is  to  use  the  techniques  of  linear  stability  theory  as  a guide  to  un  !er- 
standing  the  complex  phenomena  which  are  encountered  in  the  transition  process. 
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and  for  quantitatively  estimating  the  magnitude  and  direction  of  its  movement 
under  the  influence  of  the  various  boundary  layer  parameters.  The  effects  of 
favorable  pressure  gradient  and  surface  heating  on  boundary  layer  stability 
(and  transition)  are  automatically  included  in  classical,  paral lei -flow  stability 
theory;  however,  the  effects  of  surface  roughness  are  not  so  easily  determined. 
The  present  results  are  based  on  a phenomenological  model  for  the  effects  of 
surface  roughness  on  the  mean  flow  profiles,  and,  by  inference,  on  the  stability 
characteristics  of  a boundary  layer.  Thus,  in  the  present  work,  the  effects  of 
roughness  on  boundary  layer  transition  are  incorporated  through  a distortion  of 
the  mean  flow  profiles.  Although  the  model  is  founded  upon  experimental  infor- 
mation, that  information  is  admittedly  meager,  and  full  justification  of  the 
approach  requires  additional  experimental  verification.  The  fundamental  postu- 
late that  the  mean  flow  profiles  over  a rough  wall  are  different  from  those  which 
would  be  expected  for  a smooth  wall  is  difficult  to  refute;  the  uncertainty  lies 
in  manner  by  which  these  distortions  are  included.  Full  details  of  the  model 
are  included  in  Section  2.0;  some  representative  results  and  a qualitative 
assessment  of  the  effects  of  roughness  on  the  mean  flow  profiles  are  given  in 
Section  3.0,  and  its  effects  on  the  stability  characteristics  are  given  in 
Section  4.0. 
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2.0  APPROACH 

We  begin  by  presenting  a brief  overview  of  possible  types  of  surface  roughness 
and  then  discuss  some  experimental  results  which  show  the  mechanisms  by  which 
one  particular  type  of  roughness  affects  transition.  Building  upon  this  experi- 
mental background,  we  then  discuss  a postulated  model  for  the  effects  of  dis- 
tributed roughness  on  the  mean  flow  profiles.  Subsequent  sections  incorporate 
this  model  into  the  mean  flow,  and  finally,  present  the  results  of  calculations 
of  the  stability  characteristics  of  the  modified  mean  flow. 

2. 1 Types  of  Roughness 

Roughness  may  result  from  a variety  of  sources.  The  manufacture  of  a body  may 
require  polishing  simply  due  to  raw  materials  used;  however,  machining  itself 
may  also  be  a significant  source  of  roughness.  Additional  new  sources  of  rough- 
ness appear  once  the  body  is  introduced  into  a real  environment.  For  example, 
dust  in  air  or  comparable  particles  in  water  may  adhere  to  the  surface;  impact 
with  small  particles  may  result  in  pitting;  and  chemical  reactions  between  the 
surface  material  and  the  surrounding  fluid  may  result  in  corrosion.  Each  dif- 
ferent source  can  result  in  a different  type  of  roughness. 

In  an  attempt  to  bring  some  order  into  this  situation,  several  categories  of 
roughness  can  be  identified  for  both  analytical  and  experimental  purposes. 

These  include  a single  isolated  protrusion,  a few  isolated  protrusions  (typically 
similar  or  identical  in  structure),  and  a closely  packed  regular  or  irregular 
array  of  protrusions  (distributed  roughness).  In  addition  to  these  three- 
dimensional  types  of  roughness,  categories  of  a two-dimensional  nature,  including 
isolated  cylindrical  roughness  elements  (e.g.,  trip  wires),  and  multiple  cylin- 
drical elements  or  wavy  walls  may  also  be  encountered. 

For  each  type  of  roughness,  it  is  necessary  to  identify  those  parameters  which 
suffice  to  characterize  the  roughness.  Attempts  to  discriminate  among  roughness 
elements  of  different  shapes  have  been  reported  by  Schlichting  (1968)  for  the 
case  of  turbulent  flows  through  rough  pipes.  His  categories  include  elements 
which  are  spheres,  segments  of  sphere,  cones,  and  "short  angles."  For  the  case 
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of  distributed  roughness,  similar  categorization  of  the  shapes  of  the  roughness 
elements  is  desirable,  including  whether  the  individual  elements  are  similar  in 
shape,  or  whether  they  are  composed  of  a random  assortment  of  shapes.  However, 
some  additional  parameters  also  enter  in  the  distributed  roughness  case.  For 
example,  the  number  of  roughness  elements  per  unit  of  surface  area  is  very  impor- 
tant, as  is  the  height  distribution  of  the  elements.  Either  of  these  parameters 
could  be  used  to  define  a continuous  parametric  space  between  distributed  rough- 
ness and  single  isolated  three-dimensional  elements.  Thus,  if  the  number  density 
of  roughness  elements  were  continuously  decreased,  the  distributed  roughness 
would  approach  a single  (or,  at  least,  multiple  isolated)  roughness  element(s). 
Similarly,  if  the  height  distribution  is  continuously  changed  from  one  in  which 
all  elements  are  of  identical  heights  to  one  where  a few  elements  are  noticeably 
bigger  than  the  other  smaller  ones,  these  few  tall  elements  could  again  act  as 
isolated  roughness  elements. 

As  noted  above,  the  analysis  in  this  report  is  restricted  to  distributed  surface 
roughness,  but  in  an  attempt  to  retain  only  the  most  dominant  effects  of  the 
roughness,  the  present  analysis  characterizes  the  roughness  by  a single  param- 
eter, the  effective  roughness  height,  k.  Certainly,  this  is  an  oversimplifica- 
tion; the  number  density,  the  individual  shapes,  and  the  height  distribution  of 
the  elements  are  also  important,  but  these  details  are  ignored  for  the  present. 
Using  only  a single  parameter  to  characterize  the  roughness  should  still  give 
meaningful  qualitative  results  for  many  distributed  roughness  cases.  The  radical 
deviations  from  distributed  roughness  (and,  in  particular,  those  which  begin  to 
act  as  individual  roughness  elements)  can  not  be  treated  with  the  present  analysis. 
The  inclusion  of  their  effects  requires  additional  considerations. 

2 . 2 The  Effects  of  Surface  Roughness:  Experimental  Background 

Most  experimental  work  relating  to  the  effects  of  roughness  on  transition  has 
concentrated  on  the  relationship  between  the  various  types  of  roughness  and  the 
transition  Reynolds  number,  rather  than  on  the  microscopic  physics  in  the  neigh- 
borhood of  a roughness  element  and  how  it  affects  the  transition  mechanisms.  A 
concise  summary  of  such  results  may  be  found  in  the  review  article  by  Tani  (1969). 
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It  should  be  noted  that  even  these  "macroscopic"  experiments  are  generally  con- 
cerned with  either  a single  cylindrical  (two-dimensional)  roughness  element,  or 
an  isolated  protuberance  (three-dimensional  element),  rather  than  with  distrib- 
uted surface  roughness.  For  the  case  of  distributed  roughness,  even  the  macro- 
scopic information  concerning  the  effects  of  roughness  on  the  transition  location 
is  limited,  and  data  showing  the  effects  of  heat  transfer,  or  pressure  grad- 
ient, is  even  more  difficult  to  find. 

One  experiment  in  which  the  detailed  effects  of  roughness  on  the  movement  of 
transition  were  studied  is  the  one  by  Klebanoff  and  Tidstrom  (1972).  They 
observed  the  interaction  between  roughness  and  the  mean  velocity  profile  in  an 
experiment  wherein  a cylindrical,  two-dimensional  roughness  element  was  placed 
on  the  surface  of  a flat  plate  at  a stream-wise  location  at  which  the  boundary 
layer  was  still  laminar,  but  was  sufficiently  thick  to  submerge  the  cylinder 
completely.  Measurements  taken  downstream  of  the  cylinder  indicated  that  fluc- 
tuations inside  the  boundary  layer  were  amplified  or  damped  as  they  were  swept 
downstream,  depending  on  their  frequency.  Comparisons  with  computations  from 
flat  plate  linear  stability  theory  showed  that  the  experimentally  measured  ampli- 
fication rates  were  much  higher  than  the  values  which  would  be  predicted  from 
smooth-wall  theory;  however,  when  the  predicted  amplification  rates  were  computed 
from  the  measured  mean  velocity  profiles  (which  were  distorted  from  the  Blasius 
shape  by  the  presence  of  the  wire),  the  results  were  in  good  agreement  with  the 
measured  growth  rate. 

In  addition  to  this  increased  amplification  rate,  which  indirectly  increased  the 
fluctuation  level,  Klebanoff  and  Tidstrom  found  that  the  roughness  element  also 
served  as  a direct  source  of  disturbances  by  introducing  additional  fluctuations 
in  the  laminar  boundary  layer.  Consequently,  their  findings  can  be  summarized 
by  noting  that  roughness  affects  transition  in  the  following  two  ways: 

1.  The  presence  of  surface  roughness  alters  the  mean  velocity 
profiles  in  such  a manner  that  disturbances  in  the  laminar 
boundary  layer  are  amplified  at  a faster  rate. 

2.  The  presence  of  surface  roughness  generates  additional 
disturbances  in  the  boundary  layer  and,  hence,  changes  the 
initial  disturbance  level  before  amplification  begins. 
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These  experimental  observations  have  been  used  as  the  basis  for  our  phenomenologi 
cal  model  of  the  effects  of  distributed  roughness  on  transition  as  outlined  below 

2.3  Modeling  the  Effects  of  Distributed  Surface  Roughness 

The  experimental  evidence  of  Klebanoff  and  Tidstrom,  which  was  described  above, 
clearly  identifies  the  mechanisms  by  which  a single  two-dimensional  roughness 
element  causes  transition  to  move  forward.  Similar  detailed  experimental  infor- 
mation concerning  the  mechanisms  through  which  distributed  roughness  affects  the 
location  of  transition  is  completely  lacking,  but  it  is  to  be  expected  that  dis- 
tributed roughness  would  likewise  affect  the  location  of  transition  through  the 
same  two  basic  mechanisms:  distortion  of  the  mean  flow  profile,  and  generation 

of  additional  disturbances.  Following  the  previous  work  of  Merkle,  Kubota,  and 
Ko  (1974),  we  have  represented  the  assumed  distortion  of  the  mean  flow  by  dis- 
tributed roughness  in  terms  of  a simple  phenomenological  model.  In  the  model, 
it  is  assumed  that  distributed  roughness  elements  are  spaced  closely,  compared 
to  the  characteristic  wavelength  of  the  boundary- layer  disturbances.  Further, 
the  model  views  the  flow  over  the  many  roughness  elements  as  being  unsteady  in 
nature,  similar  to  that  commonly  observed  behind  isolated  bodies  at  intermediate 
Reynolds  numbers.  This  postulated  unsteadiness  could  occur  in  the  form  of  either 
a vortex  street  or  wake  turbulence.  In  an  integrated  sense,  these  unsteady  veloc 
ity  fluctuations  serve  near  the  wall  as  a source  of  augmented  momentum  and  heat 
transfer.  These  augmented  transfer  rates  are  modeled  by  local  eddy  di ffusi vi ties 
In  this  fashion,  the  intractable  flow  field  which  is  generated  by  the  extremely 
complicated  physical  surface  is  circumvented  by  considering  a statistical  average 
which  yields  an  "equivalent"  behavior  over  a smooth  surface.  Effectively,  a 
"turbulent  roughness  layer"  with  a thickness  of  the  effective  roughness  height 
is  assumed  to  be  imbedded  within  the  ordinary  laminar  boundary  layer. 

Quite  naturally,  this  simple  model  does  not  faithfully  reproduce  the  detailed 
flow  field  in  this  turbulent  roughness  layer;  the  complexities  of  the  flow  near 
an  actual  rough  wall  are  far  too  great  to  attempt  to  describe  completely.  The 
model  does,  however,  strive  to  incorporate  the  dominant  effects  of  the  most 
important  processes  which  occur  in  the  wall  region,  and  it  is  felt  that  these 
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will  have  a greater  influence  on  the  global  behavior  of  transition  than  the 
inconsistencies  in  the  details  of  the  flow  in  the  immediate  vicinity  of  the 
roughness.  The  model  was  originally  empirically  adjusted  to  match  the  experi- 
mental results  of  Feindt  (1957)  for  the  effects  of  distributed  roughness  on 
transition  in  the  presence  of  favorable  and  unfavorable  pressure  gradients 
(Merkle,  Kubota,  and  Ko,  1974),  and  it  has  since  been  used  to  predict  the  cor- 
rect qualitative  behavior  of  the  (individual  and  simultaneous)  effects  of  strong 
favorable  pressure  gradients,  wall  cooling,  surface  mass  addition,  and  distributed 
surface  roughness  in  compressible  boundary  layers.  Complete  justification  and 
continued  improvements  of  the  model  must  await  detailed  experimental  measurements 
of  the  flow  field  in  the  presence  of  distributed  roughness. 


The  enhanced  momentum  transfer  which  is  postulated  to  occur  near  the  surface  is 

represented  by  a momentum  diffusivity,  e^,  while  the  analogous  heat  transfer  is 

given  by  a thermal  diffusivity,  c^.  The  momentum  diffusivity  is  expected  to  be 

significant  near  the  wall,  but  it  is  expected  to  vanish  for  distances  which  are 

large  compared  to  the  roughness  height,  k.  In  mathematical  form,  the  diffusivity 

is  expressed  in  terms  of  a representative  magnitude,  e , and  a function,  F(y/k), 

max 

which  is  of  the  order  of  unity  near  the  wall,  and  goes  to  zero  for  y >>  k.  One 
such  function  which  has  this  behavior  is  the  Gaussian  function,  and  we  have 
defined  the  momentum  diffusivity  in  terms  of  it  as 


-Bjty/k)1 


= e, 


max 


Bj  =0(1] 


where  the  constant  is  required  to  be  of  order  unity  so  that  the  width  of  the 
region  of  amplified  momentum  transfer  is  similar  to  the  roughness  height.  (In 
fact,  8,  = 1.0  was  used  for  our  calculations.)  To  determine  a reasonable  esti- 
mate for  £max,  analogous  mechanisms  for  which  experimental  data  are  available 
were  considered.  Consideration  of  a turbulent  wake  [Merkle,  Kubota,  and  Ko  (1974) 
and  Schlichting  (1968)]  have  led  to 


E max 


vk  Rek 


where  the  roughness  Reynolds  number,  Re^,  is  defined  in  terms  of  the  velocity, 
u^,  and  the  kinematic  viscosity,  v^,  evaluated  at  the  roughness  height  y = k: 
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The  empirical  constant  K.'  is  expected  to  have  the  magnitude  k'  = 0[ 0.1). 
Furthermore,  it  is  expected  that  momentum  transfer  by  the  unsteady,  roughness- 
induced  velocities  has  some  threshold  Reynolds  number  below  which  the  flow  is 
completely  stable  so  that  no  "turbulent"  momentum  transfer  takes  place  in  the 
region  near  the  wall  (experimental  results  show  that  the  flow  field  over  an 
isolated  body  is  completely  stable  up  to  some  particular  Reynolds  number,  after 
which  unsteadiness  begins).  This  threshold  Reynolds  number  has  been  introduced 
into  the  model  by  allowing  the  coefficient  k'  to  have  a Van  Driest-type  dependence 
on  the  local  Reynolds  number: 

K'  = K J^l  - exp(-Rek/A+) 

Based  on  experimental  results  taken  from  bodies  in  undisturbed  flow,  it  is 
expected  that  this  threshold  Reynolds  number  is  of  the  order  of  40  (corresponding 
to  the  Reynolds  number  at  which  the  vortex  street  behind  a cylinder  begins). 


To  summarize,  the  momentum  diffusivity  due  to  roughness  is  given  by 


l'M 


K v.  Re 
e k 


I 

k ) 


1 - 


-RVA+t 


-3j(y/k)? 

e 


For  all  the  calculations  reported  here,  the  three  constants  in  the  model  were  set 
equal  to  their  order  of  magnitude  estimates  described  above.  Thus,  these  con- 
stants were  given  the  values. 


K = .094 

e 

A+  = 40 
61  = 1 


Finally,  the  thermal  diffusivity  due  to  the  surface  is  defined  in  terms  of  a 
turbulent  Prandtl  number,  P^, 


f H = £ M/Pf 
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For  all  calculations  reported  herein,  the  turbulent  Prandtl  number  was  given 
the  value  of  unity.  The  values  for  these  four  constants  are  identical  to  those 
chosen  by  Merkle,  Kubota,  and  Ko  (1974)  and  by  Merkle  (1976)  for  compressible 
rough  wall  boundary  layers.  None  of  the  constants  nor  any  of  the  assumptions  in 
the  model  have  been  changed  to  apply  the  analysis  to  boundary  layers  in  water. 

The  total  or  effective  viscosity  and  thermal  conductivity  are  given,  respectively, 
by 

PT  = P + peM 

kT  = k + PCpeH 

Here  p and  k are  the  usual  viscosity  and  thermal  conductivity,  p is  the  density, 
and  Cp  is  the  specific  heat.  The  effect  of  the  enhanced  diffusivities  is  to 
decrease  the  velocity  near  the  wall  and  to  make  the  temperature  profile  thicker. 
Qualitative  graphs  illustrating  this  are  shown  in  Figure  1.  The  total  diffu- 
sivities, reflecting  the  sum  of  the  laminar  diffusivity  and  an  "enhanced"  contri- 
bution due  to  the  presence  of  surface  roughness,  are  shown  schematically  in 
Figure  2.  Near  the  wall  the  total  diffusivity  is  large,  but  at  heights  greater 
than  the  roughness  height  the  usual  laminar  value  is  regained.  Note  from  Figure  1 
that  when  the  velocity  near  the  wall  decreases  sufficiently,  the  mean  profile 
develops  a point  of  inflection  and  is  inviscidly  unstable. 

The  above  equations  and  diffusivities  provide  all  the  additional  information 
needed  to  incorporate  the  effects  of  distributed  surface  roughness  into  the  cal- 
culation of  the  distorted  mean  profiles  in  water,  including  the  presence  of 
heating  and  pressure  gradients.  Calculations  of  this  nature  are  described  in 
the  following  section,  and  the  subsequent  section  contains  the  results  of  a 
stability  analysis  of  the  deformed  profiles. 
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3.0  MEAN  FLOW  ANALYSIS 


As  noted  by  Klebanoff  and  Tidstrom  (1972),  one  of  the  mechanisms  by  which  rough- 
ness influences  transition  is  through  the  distortion  of  the  mean  flow  profiles. 
Having  discussed  how  the  effects  of  distributed  surface  roughness  are  to  be 
modeled,  we  proceed  to  the  development  of  the  specific  equations  used  for  the 
mean  flow  and  indicate  the  results  attained. 

3. 1 Theory 

The  equations  for  a compressible,  laminar  boundary  layer  about  an  axisymmetric 

body  of  radius  r are: 
o 

Conservation  of  mass 
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Note  that  the  usual  viscosity,  p,  and  conductivity,  k,  have  been  replaced  by 
their  rough-wall  counterparts,  Uy  and  ky.  The  temperature,  T,  which  is  needed 
to  evaluate  the  viscosity  and  thermal  conductivity,  can  be  calculated  from  the 
stagnation  enthalpy,  H,  according  to 

H = C (T  - T ,)  + ‘2u2 
p ' ref' 

since  the  specific  heat  of  water  is  essentially  constant  in  the  region  of  inter- 
est.  (Trep  is  merely  a reference  temperature  which  does  not  affect  the  calcula- 
tions. ) 

Although  these  partial  differential  equations  could  be  solved  directly  (using  an 
appropriate  numerical  procedure),  their  solutions  would  depend  on  the  specific 
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body  shapes  and  pressure  distributions  which  were  used,  whereas,  for  these 
initial  results,  we  are  more  interested  in  obtaining  a qualitative  description 
of  the  movement  of  the  transition  location  in  response  to  various  heights  of 
surface  roughness.  Accordingly,  we  have  chosen  to  work  with  the  similarity 
equations  rather  than  the  full  partial  differential  equations.  The  similar 
form  of  the  equations  are  somewhat  simpler  to  solve  (since  they  are  ordinary 
differential  equations),  but  they  also  make  it  much  easier  to  parameterize  the 
numerous  variables  in  the  problem,  and,  thus,  to  obtain  a more  rapid  understanding 
of  the  global  effects  of  distributed  roughness  on  transition,  as  predicted  by 
the  model.  Future  calculations  are  planned,  however,  in  which  the  onset  of 
transition  on  realistic  bodies,  will  be  predicted.  For  these  future  calculations, 
the  full  partial  differential  equations  will  be  solved  by  means  of  a standard 
numerical  method. 

It  is  worth  noting  that  the  axisymmetric  nature  of  the  geometry  of  the  underwater 
body  enters  the  analysis  automatically  through  the  mean  flow  equations.  This 
results  in  significant  differences  between  axisymmetric  and  planar,  two- 
dimensional  boundary  layers  (e.g.,  the  growth  rate  of  the  boundary- layer  thick- 
ness). By  contrast,  it  is  an  established  fact  that  the  form  of  the  stability 
or  disturbance  equations  (discussed  in  the  next  section)  is  the  same  for  axi- 
symmetric  and  planar,  two-dimensional  flows;  however,  the  coefficients  differ 
since  they  are  derived  from  different  mean  flow  equations. 

To  determine  the  equations  for  similarity  solutions  corresponding  to  axisymmetric 
flow,  we  make  use  of  the  Levy-Lees  transformation  (Hayes  and  Probstein,  1955) 

£(x)  = f pQpQuQr^dx 
Jo  e e e o 


n(x,y) 


Vo_  n 
yfe  Jo 


P dy 


The  subscript  "e"  here  is  used  to  denote  evaluation 
layer.  For  axisymmetric  flow  there  exists  a stream 
similarity  solution,  can  be  written  as 


at  the  edge  of  the  boundary 
function,  H',  which,  for  the 
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Using  the  standard  methods,  the  velocity  can  then  be  computed  from 

u = ugF ' (n), 

where,  here  the  prime  denotes  differentiation  with  respect  to  eta.  Similarly, 
the  total  temperature  can  be  expressed  in  terms  of  the  similarity  function,  G, 
as 

G (n)  = H/He. 

Then,  by  introducing  the  local  non-dimensional  viscosity  and  the  effective 
Prandtl  number, 

C = Py/Pg*  Py  ” CpPy/ky 

along  with  the  Falkner-Skan  pressure  gradient  parameter, 

du_ 
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x _?L 
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the  relevant  equations  become 
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We  are  able  to  make  several  simpl ications  based  on  the  fact  that  we  are  interested 
in  low-speed  flow  in  a heated  water  boundary  layer: 


water  is  essentially  imcompressible  and  p 


dH« 


• the  free  stream  enthalpy  is  constant  and  hence  . - = 0 

• because  u is  small  and  the  specific  heat  of  water  is  large, 
Ug/He  <<  1 for  only  moderate  heating. 

The  system  of  equations  actually  used  to  generate  the  mean  flow  profiles  thus 
simplifies  to 
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, 1 

(CF1  1 )'  + FF"  + e[l  - (F ' ) 2]  = 0 
G'  )’  + FG'  =0 

The  first  of  these  equations  is  coupled  to  the  second  through  the  dependence  of 
viscosity  and  the  roughness  model  on  temperature.  The  system  is  completed  by 
specifying  the  boundary  conditions  at  the  wall, 

no-suction,  no-slip  conditions:  F ( 0 ) = F 1 ( 0 ) = 0 

prescribed  wall  temperature:  G(0)  = G ; 

w 

and  the  free  stream  values  which  are  approached  at  the  edge  of  the  boundary  layer, 

F'  (n)  ->■  1,  G(n)  -*■  1 as  ri  + “■ 

These  equations  have  been  integrated  numerically  by  means  of  a standard  Runge- 
Kutta  technique. 

We  note  in  passing  that  these  boundary  conditions  for  the  rough-wall  problem 
are  identical  to  those  for  the  smooth  wall  problem  (i.e.,  they  are  specified  at 
y = 0 rather  than  on  the  actual,  irregular  surface).  Since  we  are  working  with 
boundary  layer  equations,  such  an  approximation  at  first  glance  appears  accept- 
able; it  is  acceptable  to  treat  an  irregular  surface  as  a plane,  so  long  as  the 
radius  of  curvature  is  larger  compared  to  the  boundary  layer  thickness.  However, 
closer  inspection  reveals  that  the  local  radius  of  curvature  of  the  surface  is  of 
the  order  of  the  roughness  height,  so  that  not  only  the  boundary  layer  curvature 
corrections,  but  also  the  complete  elliptic  equations,  should  be  used  in  this 
region.  The  requirement  for  the  elliptic  equations  is,  of  course,  a statement 
that  the  pressure  is  no  longer  impressed  by  the  outer  flow  field,  but  that  it 
becomes  a local  variable  which  must  be  allowed  to  vary  both  along  and  across  the 
boundary  layer.  Despite  these  many  small  scale  fluctuations,  the  spirit  of  the 
present  analysis  is  to  attempt  to  include  the  global  effects  of  these  local  pres- 
sure perturbations  through  a phenomenological  model  in  much  the  same  way  that 
local  (turbulent)  pressure  fluctuations  in  a turbulent  boundary  layer  are  ignored 
in  the  computation  of  the  mean  flow. 

For  completeness,  we  note  that  the  physical  properties  of  water  which  were  used 
in  the  present  analysis  were  computed  from  the  formulas  suggested  by  Lowell  and 
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Reshotko  (1974).  After  a review  of  the  literature  for  empirical  expressions  for 
the  dependence  of  the  viscosity  and  thermal  conductivity  of  pure  water  on  temper- 
ature, they  recommended  the  following  expressions: 


(h/u2ooc^ 


-2.302585 
T - 164 


1.37023  + 8.36  x 10"4  (T 


293) 


(T  - 293) 


and 

k = 104  <j-9. 90109  + 0.1001982  T - 1.873892  x 10'4  T2  + 1.03957  x 10'7  T3j\ 

where  the  reference  viscosity  is  h2u°C  = These  expressions  have  been 

chosen  so  that  the  first  derivatives,  as  well  as  the  functions  themselves,  can  be 
represented  with  a high  degree  of  accuracy.  In  these  formulas,  the  temperature 
T is  prescribed  in  degrees  Kelvin,  the  molecular  viscosity  y is  given  in  centipoise, 
and  the  thermal  conductivity  k is  given  in  ergs  per  centimeter-second-degree 
Kelvin. 


The  effects  of  roughness  enter  through  the  use  of  Vj  and  kT  in  place  of  y and  k, 
respectively.  By  setting  the  roughness  height  to  zero  (or  to  a very  small  value) 
and  equating  the  free  stream  and  wall  temperatures,  it  was  possible  to  generate 
the  conventional  (adiabatic)  Falkner-Skan  solutions  as  a verification  of  the 
reliability  of  the  code  used  to  solve  the  similarity  equations.  In  order  to 
verify  the  code  when  there  is  a heat  flux  across  the  boundary  layer,  use  was  made 
of  a relation  between  G and  F for  a (smooth)  flat  plate  (i.e.,  6=0)  established 
by  Pohlhausen  (Curie,  1962).  If  eg)  = ( T - Tw)  / (T  - T ),  where  Tw  is  the  wall 
temperature,  then  ® is  given  in  terms  of  the  Blasius  profile  F and  the  free  stream 
Prandtl  number  Pr  by 


® (n) 


I 


Pr 


dn 


The  function  aQ(Pr),  which  has  been  tabulated,  is  given  by 

a0(pr)  = |/  exp|-Pr  J Fdq)  dn  | 

and  is  reasonable  well  approximated  by 

aQ(Pr)  = 0.664(Pr)1/3. 
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3.2 


Results 


A number  of  mean  velocity  and  temperature  profiles  were  calculated  for  parametric 
values  of  heating,  pressure  gradient,  wall  temperature,  and  effective  roughness 
height.  The  overall  results  found  are  exemplified  by  the  profiles  for  a flat 
plate  (0  = 0)  which  appear  in  Figure  3.  The  boundary  layer  thickness  6 was  taken 
to  be  the  value  of  y where  u/uQ  = .999.  Normalized  velocities  u/u„  and  temper- 

G 0 

atures  ® = AT/AT  are  shown  as  function  of  y/6  for  various  values  of  the  non- 
w 

dimensional  roughness  height,  k/0.  Here  AT=T-T,AT  = T - T , k is  the 
3 e w w e 

effective  roughness  height,  and  0 is  the  momentum  thickness  of  the  distorted  mean 
profile.  The  results  quantitatively  verify  the  qualitative  profiles  anticipated 
in  Figure  1. 

As  roughness  increases,  the  gradients  of  both  the  velocity  and  temperature  pro- 
files decrease  near  the  wall  as  a consequence  of  the  increased  diffusivity  there. 
Near  the  boundary  layer  edge,  the  profiles  are  not  influenced  greatly  by  the 
diffusi vities.  This  tends  to  introduce  an  inflection  point  into  the  mean  velocity 
profile,  making  the  profile  more  unstable  from  an  inviscid  point  of  view.  Because 
of  the  large  Prandlt  number  of  water  [Pr  =$(10)],  the  smooth  wall  thermal  bound- 
ary layer  is  considerably  thinner  than  the  velocity  boundary  layer.  Since  the 
"turbulent  Prandlt  number"  P^.  is  of  order  unity,  the  effect  of  roughness  is 
experienced  equally  by  both  the  velocity  and  temperature  profiles  near  the  wall. 
The  consequence  of  this  is  that  a moderate  degree  of  roughness  tends  to  drive 
the  thickness  of  both  the  velocity  boundary  layer  and  the  temperature  boundary 
layer  towards  each  other. 

Note  that  for  the  flat  plate  case,  the  k/0  = 0 and  k/0  = .3  the  curves  are  essen- 
tially coincident.  This  suggests  that  a certain  range  of  roughness  may  be  con- 
sidered equivalent  to  a smooth  wall.  For  k/0  = 1 the  curves  deviate  modestly 
from  the  smooth  wall  profiles,  and  for  k/0  = 2 they  deviate  substantially.  How- 
ever, stability  calculations  presented  in  the  next  section  demonstrate  that  only 
a modest  deviation  from  the  smooth  wall  profile  is  needed  to  result  in  a drastic 
change  in  the  stability  characteristics. 
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4.0  STABILITY  ANALYSIS 


In  the  previous  section  we  presented  a model  to  evaluate  how  distributed  surface 
roughness  alters  the  mean  flow  profiles.  These  distorted  profiles  should  have 
different  stability  characteristics  which  will  influence  the  onset  of  the  transi- 
tion phemomena.  The  transition  criteria  which  we  plan  to  utilize  in  future  analy- 
ses, described  below,  is  primarily  based  upon  spatial  amplification  rates. 

It  is  appropriate  to  recall  here  an  alternative  transition  mechanism,  other  than 
the  deformation  of  the  mean  profiles,  as  it  relates  to  stability  theory.  As 
previously  noted,  Klebanoff  and  Tidstrom  noted  that  roughness  may  introduce  addi- 
tional disturbances  into  the  flow,  changing  the  spectral  content.  If  we  postulate 
that  transition  occurs  when  a critical  amplitude  has  been  attained,  it  follows 
that  when  the  initial  amplitudes  present  are  larger,  transition  will  occur  for  a 
smaller  value  of  the  amplification  ratio.  This  could  be  incorporated  into  "e11" 
transition  criterion  by  connecting  the  appropriate  choice  of  "n"  to  the  nature  of 
the  disturbance  spectrum.  We  believe  that  this  does  not  violate  the  spirit  in 
which  Smith  and  Gamberoni  (1956)  intended  the  criteria  to  be  used.  The  calcula- 
tions involved  in  an  "e11"  type  prediction  derive  directly  from  stability  theory. 

4. 1 Theory  and  Procedure 

The  flow  in  a boundary  layer  can  be  decomposed  into  the  sum  of  a mean  flow  and 
disturbances.  If  the  disturbances  are  sufficiently  small,  the  relevant  equations 
can  be  linearized.  If,  further,  one  makes  a parallel  flow  assumption,  i.e.,  that 
the  mean  profiles  and  the  pressure  gradient  vary  slowly  in  the  axial  direction, 
then  the  disturbances  can  be  decomposed  into  components  of  the  form 

<J>'  (x,y,t)  = $(y)ei(ax-wt) 

(Recall  that,  from  a structural  point  of  view,  the  disturbance  equations  for  axi- 
symmetric  and  planar,  two-dimensional  flows  are  identical,  differing  only  in  the 
evaluation  of  the  coefficients  from  the  mean  flow  profiles.)  This  results  in  a 
system  of  ordinary  differential  equations  for  functions  of  y.  To  study  spatial 
stability,  which  is  the  focus  of  this  report,  the  frequency  w is  taken  to  be  real 
and  prescribed.  Since  the  disturbance  or  stability  equations,  as  well  as  the 
boundary  conditions,  are  linear  and  homogeneous,  nontrivial  solutions  are  only 
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possible  for  complex  wave  numbers  a which  are  eigenvalues  of  the  system.  Writing 

a = + i a.j 

we  can  express  the  local  spatial  amplification  rate  for  a disturbance  of  a given 
frequency  as  (-a,).  Thus,  once  the  mean  profiles  are  known,  it  is  possible  to 
integrate  the  linear  stability  equations  at  each  stream-wise  location  (or  Reynolds 
number)  to  obtain  the  local  amplification  rate  associated  with  each  individual 
frequency  component. 

The  indicated  integration  has  been  carried  out  using  a numerical  code  developed 
by  Lowell  and  Reshotko  (1974),  after  making  appropriate  modifications  so  that  it 
accepts  the  mean  flow  profiles  we  have  generated.  Modifications  to  include  terms 
introduced  into  the  stability  equations  by  the  enhanced  viscosity  have  not  been 
made  at  the  present.  It  is  anticipated  that  their  effects  on  the  results  would 
be  smaller  than  those  caused  by  the  distortions  of  the  mean  flow  profile.  The 
code  uses  an  iterative  method  and  incorporates  an  orthogonal ization  procedure  to 
minimize  parasitic  error  growth.  The  equations  include  temperature  fluctuations 
and  reflect  the  variation  of  viscosity  and  thermal  conductivity  with  temperature, 
using  the  analytic  expressions  indicated  in  the  previous  section. 

Nondimensional  frequencies  are  prescribed: 

v 

n - e 

SZ  = CO  5 — 

ue 

Eigenvalues  a were  calculated  as  a function  of  fi,  Re0,  0,  ATw,  and  k/0  (non- 
dimensional  roughness  height).  This  constitutes  a five-dimensional  parameter 
space,  and  so  it  was  necessary  to  choose  values  carefully  in  order  that  results 
could  be  attained  within  reasonable  time  and  financial  constraints.  For  quali- 
tative assessment,  an  ambient  water  temperature  of  59°F  was  used  for  all  calcu- 
lations. Values  of  6 chosen  were  -0.05  (adverse  pressure  gradient,  dp/dx  > 0), 

0.0  (flat  plate,  dp/dx  = 0),  and  0.2  (favorable  pressure  gradient,  dp/dx  < 0). 

Rer  was  taken  as  an  indicator  of  position  and  was  held  constant  for  each  0.  To 
illustrate  the  relative  insensitivity  resulting  from  this  choice  over  Rer*  or  Re . , 

6 O 

plots  of  H = 6*/0  and  6/0  are  given  for  the  flat  plate  as  a function  of  k/0  in 
Figure  4.  Values  of  Re0  were  chosen  so  that,  for  the  adiabatic  smooth  wall.  Re 
(or  Re.*)  would  lie  approximately  in  the  region  of  maximum  instability.  [(Ref- 
erence can  be  made  to  the  spatial  stability  maps  generated  by  Wazzan,  Okamura, 
and  Smith  (1968).] 
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The  values  are 


3 

Refl 

Re6* 

-0.05 

224 

600 

0 

385 

1000 

0.2 

2103 

5071 

0,  k/0  = 0) 


As  noted  in  the  introduction,  the  calculations  performed  are  intended  to  be  used 
for  transition  prediction.  A cumulative  amplification  ratio  can  be  calculated 
from  the  stability  analysis  according  to 


A/Aq  = exp  j 


V 


where  Ao(w)  and  A(w)  are,  respectively,  the  initial  and  local  amplitudes  of  a 
disturbance  with  frequency  to.  The  en  transition  criterion  which  has  been  developed 
by  Smith  and  Gamberoni  (1956)  predicts  that  the  onset  of  transition  will  occur  when 
this  amplification  ratio  has  reached  a critical  level 


for  some  frequency.  Previous  experience  with  adiabatic  boundary  layers  indicates 
that  the  parameter  n is  generally  of  order  9 or  10. 


4.2  Results 

The  principal  results  of  the  stability  calculations  are  illustrated  in  Figures  5 
through  10,  which  show  the  dependence  of  the  nondimensional  amplification  rates 

(-a.0)  as  a function  of  ft,  the  frequency,  for  3=0  (Re„  = 385),  6 = 0.2  (Re  = 

l o 

2103),  and  8 = -0.05  (Ret>  = 224).  The  figures  occur  in  pairs  for  each  value  of  8. 
The  first  shows  curves  for  a variety  of  nondimensional  roughness  heights  when 
ATw  = 0°F;  the  second  of  each  pair  shows  the  same  for  both  AT  = 0°F  (solid  lines) 
and  AT  = 30°F  (dashed  lines). 

Insofar  as  the  curves  are  qualitatively  similar  in  shape,  it  seem  reasonable  to 
characterize  each  by  the  peak  amplification  rate  (“aj°)max-  With  this  simpli- 
fication, there  still  remain  three  parameters  whose  interrelations  we  wish  to 
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study,  viz.  k,  ATw,  and  3 (or  dp/dx).  We  accomplish  this  by  fixing  3 and  then 
generating  two  types  of  plots  for  each  3:  curves  showing  (-a. 6)  as  a function 

of  k/6  for  various  values  of  ATw  (Figures  11,  12,  and  13),  and  curves  showing 
(-a-0)  as  a function  of  AT  for  fixed  values  of  k/0  (Figures  15,  16,  and  17). 

Figures  5,  7 and  9 demonstrate  the  anticipated  result  that  the  level  of  instabil- 
ity increases  as  the  roughness  increases  in  the  absence  of  heating.  It  is  also 
worthwhile  noting  that  the  range  of  unstable  frequencies  expands  as  roughness  is 
increased.  When  heating  is  present  (Figure  6,  8,  and  10)  a smooth  wall  is  stabi- 
lized to  some  degree,  as  previously  known.  However,  roughness  plays  a significant 
role  in  modifying  this  behavior.  Initially  there  is  a range  of  values  for  which 
a rough  wall  is  indistinguishable  from  a smooth  wall;  a further  increase  in  rough- 
ness then  results  in  a moderate  lessening  of  the  stablizing  effects  of  added  heat. 
However,  once  the  roughness  height  surpasses  a critical  value,  heating  exhibits 
a marked  destabilizing  effect  on  the  boundary  layer.  Figures  11,  12,  and  13 
present  these  results  in  a more  succinct  fashion.  From  these  figures  it  is  pos- 
sible to  interpolate  values  of  k/0  which  correspond,  qualitatively,  to  negligible, 
moderate,  and  drastic  changes  in  the  stability  characteristics.  Appropriate 
values  are  tabulated  in  Figure  14.  These  data  permit  us  to  draw  some  conclusions 
about  the  influence  of  a pressure  gradient  in  the  presence  of  roughness.  Specifi- 
cally, the  more  favorable  the  pressure  gradient  (i.e.,  as  3 increases),  the 
narrower  is  the  range  of  roughness  heights  for  which  a wall  may  be  treated  as 
smooth.  Also,  the  more  favorable  the  pressure  gradient,  the  smaller  are  the  cri- 
teria for  what  constitutes  "critical  roughness"  beyond  which  there  is  a drastic 
tendency  to  destabilize. 

Figures  15,  16,  and  17  indicate  that,  for  a given  pressure  gradient  and  roughness 
height,  there  exists  an  "optimal"  amount  of  heating  which  stabilizes  the  boundary 
layer  as  much  as  possible,  and  a maximum  amount  beyond  which  heating  only  destabi- 
lizes. Clearly  the  former  parameter  is  of  principal  interest.  For  sufficiently 
large  roughness,  the  "optimal"  amount  of  heating  turns  out  to  be  no  heating  at 
all. 
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5.0  SUMMARY  AND  CONCLUSIONS 


An  analysis  has  been  conducted  to  estimate  the  interrelations  among  the  effects 
of  heating,  shaping,  and  particularly  distributed  surface  roughness  on  the  tran- 
sition characteristics  of  high  Reynolds  number  boundary  layers  in  water.  The 
estimates  of  the  behavior  of  the  transition  location  have  been  based  upon  the 
results  of  linear  stability  theory.  The  effects  of  distributed  surface  rough- 
ness have,  in  turn,  been  included  in  the  stability  analysis  by  means  of  a phenom- 
enological model  for  the  effects  of  roughness  on  transition.  The  model  visualizes 
that  roughness  affects  the  location  of  transition  by  causing  a distortion  in  the 
mean  flow  profile  (from  its  smooth-wall  shape)  in  such  a manner  as  to  alter  its 
stability  characteristics.  . 

The  model  which  has  been  used  for  these  incompressible,  water  boundary  layer 
predictions  is  identical  to  one  previously  developed  and  used  for  compressible, 
air  boundary  layers.  For  the  case  of  compressible  boundary  layers,  the  model 
has  been  demonstrated  to  predict  the  correct  qualitative  trends  for  the  effects 
of  rougness  on  transition,  as  well  as  for  the  individual  and  simultaneous  effects 
of  streamwise  pressure  gradients,  surface  heat  transfer,  surface  ablation,  and 
Mach  number,  in  the  presence  of  roughness.  Similar  predictions  for  the  mutual 
interaction  between  surface  roughness,  surface  heat  transfer  and  pressure  gradient 
in  water  boundary  layers  are  given  in  the  present  report.  The  present  results 
indicate  that  the  effects  of  these  other  parameters  on  transition  can  be  quali- 
tatively different  in  the  presence  of  roughness  as  compared  to  their  behavior  on 
smooth  walls.  (As  an  example,  heating,  which  is  strongly  stablizing  in  the 
presence  of  smooth  walls,  becomes  strongly  destablizing  on  rough  walls.)  Although 
experimental  results  to  validate  these  predictions  are  presently  lacking,  it  is 
noted  that  the  model  predicts  similar  reversals  in  the  effects  of  a given  param- 
eter on  transition  in  the  presence  of  rough,  as  compared  to  smooth,  walls  for 
air  boundary  layers,  and  these  predictions  have  been  supported  by  a large  bulk 
of  experimental  observations. 

Some  specific  conclusions  and  observations  about  the  roughness  model  and  about 
the  present  predictions  of  the  effects  of  roughness  on  transition  are  outlined 
below: 
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1.  Roughness  represents  a powerful  destabilizing  influence  in  a bound- 
ary layer.  The  roughness  model  predicts  that  there  is  a range  of  roughness 
heights  for  which  boundary  layers  behave  almost  as  if  they  were  on  a smooth  wall, 
but  above  this  "threshold"  roughness  level,  the  effects  of  roughness  quickly  grow 
to  where  they  completely  dominate  the  stability  and  transition  characteristics  of 
the  boundary  layer. 

2.  In  the  presence  of  favorable  pressure  gradients,  this  threshold 
roughness  level  (when  measured  in  terms  of  k/6)  is  considerably  smaller  than  for 
a flat  plate.  This  trend  continues  to  unfavorable  pressure  gradient  boundary 
layers  where  the  threshold  roughness  level  is  greater  than  for  a flat  plate. 

3.  As  the  roughness  height  is  increased,  the  stabilizing  effects  of 
both  favorable  pressure  gradient  and  heating  on  the  growth  of  disturbances  in  the 
boundary  layer  is  diminished,  and  when  the  surface  roughness  is  sufficiently 
large,  the  effects  of  heating  no  longer  serve  to  stabilize  the  boundary  layer, 
but  actually  serve  to  destabilize  it. 

4.  The  present  model  for  the  effects  of  distributed  surface  roughness 
on  the  stability/transition  characteristics  of  a boundary  layer  must  certainly 
be  tested  against  a well-controlled  laboratory  experiment  which  clearly  identi- 
fies the  mechanisms  whereby  roughness  affects  transition.  Without  the  avail- 
ability of  such  a check  case,  there  are  a number  of  minor  improvements  which 
could  be  made  to  the  model.  For  example,  the  effects  of  the  initial  disturbances 
which  are  generated  by  the  roughness  elements  could  be  included,  and  the  effects 
of  the  enhanced  viscosity  could  be  incorporated  directly  into  the  stability  cal- 
culations. The  magnitude  of  such  effects  is,  however,  expected  to  be  minor  and 
is  not  considered  to  have  major  effects  of  the  movement  of  the  transition  loca- 
tion due  to  roughness. 

5.  Finally,  it  is  noted  that  the  present  roughness  model  is  applicable 
only  to  distributed  surface  roughness  and  cannot  be  applied  (without  modification) 
to  predict  the  effects  of  individual  three-dimensional  roughness  elements,  or 
two-dimensional  roughness  of  the  "wavy-wall"  type. 
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Figure  7.  Rough  Wail  Amplification  Rates, 
Favorable  Pressure  Gradient,  (3  = .2  Rea  = 2103 
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Figure  14.  Effect  of  Roughness  on  Amplification  Rates 
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